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2 ' Abstract: Liquid in grooved capillaries, made by e.g. inserting a plate 

O I in a cylindrical tube, exhibits unusual spreading and flow properties. One 

0\ . example is capillary rise, where a long, upward tongue on top of the usual 

^ • meniscus has been observed along the groove. We attribute the underlying 

C ' mechanism to a thermodynamic instability against spreading for a (partial 

'^ I or complete wetting) liquid in a sharp groove whose opening angle a is less 

^ ! 

O . than a critical value ac = tv — 29. The equilibrium shape of the tongue 

CJ ■ 

is determined analytically. The dynamics of liquid rising is studied in the 

Lj viscous regime. When the diameter of the tube is smaller than the capillary 

H ' length, the center part of the meniscus rises with time t following a t^' ^-law, 

while the tongue is truncated at a height which grows following a t^/"^-law. 
Sharp groove also facilitates release of gas bubbles trapped inside a capillary 
under the action of gravity. 

PACS numbers: 68.10.Gw, 82.65.-i, 87.45. Ft 



1. Introduction 

It is well-known that capillary rise in a thin tube is inversely proportional to the 
diameter of the tube[l]. This law holds for partial as well as complete wetting. In 
this paper we show that, when a sharp, vertical groove of sufficiently small opening 
angle is constructed in the tube, a liquid tongue of macroscopic thickness appears 
inside the groove. The tongue extends to arbitrarily high altitudes, independent of 
the diameter of the tube. 

Our work started with the need to fill thin glass tubes with certain kind of ionic 
solution, to be used as microelectrodes in physiological experiments. A practical 
problem is to avoid trapping of gas bubbles in the tube. When tubes with a circular 
cross-section are used [see Fig. 1(a)], it is very difficult to get rid of bubbles particu- 
larly in the sharp tip region. A simple solution was found by inserting a plate inside 
the raw tube before stretching it into the final shape using standard glass-making 
technique [see Fig. 1(b)]. Tubes with cross-sectional shapes shown in Figs. l(c)-l(e) 
were found to be equally good for this purpose[2]. 

The tubes which do not trap gas bubbles have a common feature that there are 
sharp grooves inside. When such a tube is held vertically, liquid on top of a trapped 
bubble develops tongues which extend down along these grooves. As a result, trapped 
bubbles rise easily to the top under gravity. The spontaneous flow of liquid along 
the grooves suggests that there is no free energy barrier for spreading. Indeed, when 
these tubes (about 1 mm in diameter and initially empty) are placed vertically in 
contact with the solution, we discover that, on top of the normal capillary rise, the 
liquid creeps up along the grooves to very high altitudes. 

The new physics introduced by the groove geometry is the increased tendency for 
the liquid to spread. The bottom part of the grooves which appear in the construc- 
tions shown in Figs. l(b)-(e) can be approximated by a wedge of opening angle ot. 
As we shall show below, within the classical thermodynamic theory, one can identify 
a critical angle 

«c = 7r-26', (1) 

where Q is the equilibrium contact angle. For a < olc-, spreading occurs even when 
the liquid does not wet the wall completely. Depending on the initial state of the 
liquid, gravity may either enhance or act against spreading. In this paper we focus on 
capillary rise against gravity. In this case there is an equilibrium liquid-air interface 



profile which can be determined analytically. The dynamics of the rising process will 
also be considered. 

The paper is organized as follows. In section 2 we discuss the static aspects of 
the problem. Section 3 contains an analysis on the dynamics of the rising process in 
the viscous regime. Discussion and conclusions are presented in Sec. 4. 

2. The Equilibrium Problem 

2.1. Critical angle for spreading 

In the absence of gravity, the instability of a liquid drop against spreading in 
a sharp groove with an opening angle a < Uc can be seen as follows. Consider a 
(fictitious) liquid column of length L that fills the bottom part of the groove. The 
cross-section of the column is taken to be an isosceles triangle with the two equal 
sides on the wall and a perpendicular of length d. For L ^ d, the total surface free 
energy of the liquid is given by 

F = 2rfL7tan(a/2) + 2dL{-fsL - 7sg) sec(a/2), (2) 

where 7, 750, and jsl are the surface tensions of liquid and gas, solid and gas, and 
solid and liquid, respectively. Using the Young's formula 

7COS6' = 75G -75L, (3) 

Eq. (2) may be rewritten as 

F = 2dL7[sin(a/2) - cos^]/ cos(a/2). (4) 

Equation (4) shows that, for sin(Q;/2) < cos 9 or a < ac = n — 26 and a fixed volume 
V = (i^Ltan(Q;/2), the free energy of the column decreases without bound as L ^ 00. 
We are thus led to the conclusion that any liquid drop of finite extent is unstable 
against spreading for a < ac- 

2.2. Equilibrium shape of the tongue 

When the groove is placed vertically, spreading is countered by gravity. In a 
cylindrical tube, the macroscopic rise z of the liquid is limited to a finite value even 
in the case of complete wetting. In contrast, for tubes with sharp grooves, the finite 



free energy gain per unit length in the groove [see Eq. (4)] is able to overcome the 
gravitational energy pgz per unit volume, provided the product dz is sufficiently 
small. Thus a thin tongue of the liquid can climb up to an arbitrary height along 
the bottom of the groove. The width d of the tongue, however, should be inversely 
proportional to the height. Note that the tongue is not a precursor film. It is 
completely describable within the classical thermodynamic framework. 

To determine the equilibrium profile of the liquid-air interface z{x, y), one has to 
solve the equation [1] 

pgz = 2-fH, (5) 

subjected to the boundary condition that the contact angle of the liquid at the wall is 
9. Here p is the mass density of the liquid (or more precisely, the difference between 
the mass density of the liquid and that of the gas), and H is the mean curvature of 
the interface. To obtain a complete solution to (5) for the geometries illustrated in 
Fig. 1, one has to resort to numerical means. Here we are mainly interested in the 
profile of the tongues. The problem can then be simplified by considering liquid rise 
in an infinite wedge formed by two planar vertical walls at an angle a, as shown in 
Fig. 2(a). 

Anticipating a slow variation of the cross-sectional shape of the tongue with 
the height at high altitudes, one may, in a first approximation, ignore the interface 
curvature in the vertical direction. According to Eq. (5), the horizontal cross-section 
of the tongue should then take the form of a dart with a circular base, as illustrated 
by the closed area OCBDO in Fig. 2(b). The radius R of the arc is given by 

R = a^/{2z), (6) 

where a = [2'y/(pg)]^''^ is the capillary length. The angle between the tangent of the 
circle and the wall at the point of intersection is 9. From the geometry one obtains 
the following equation relating the distance d = OB and the radius R of the circle, 

d = fl(^-l). (7) 

Combining (7) with (6), we see that the width of the tongue is inversely proportional 
to the height when the opening angle a is less than etc- 

It turns out that the leading order correction to (7) can also be calculated ana- 
lytically when the ratio R/z = \{a/zY is small. For this purpose let us introduce 



polar coordinates (r, (p) in the horizontal xy-p\a,ne centered at A, 



X = X — r cos 1 



y = r sin > 



(8) 



where X = d + R is the distance between O and A. The liquid-gas interface position 
is specified, for small R/z, by 



r(d>, z) = R 






(9) 



The mean curvature H can now be expressed in e. After some algebra the final result 
is given by 



H 



1 



l-(-Y(^ + e + l-x' + l{l-Xoos<pr 



2Rl \zJ 



+ 



(10) 



where x = X/R = cos9/ sm(a/2). Inserting (10) into (5) and using (6), we obtain, 
to leading order in (R/z)'^, 

+ e + 1 - x' + ^(1 - X cos <^)2 = 0. (11) 

Equation (11) has a general solution which respects the 4> -^ —(p symmetry, 

5 3 1 

e((/)) = C cos (p - - + -x(psm(p + -x^ cos^ (p. (12) 



The constant C in (12) is to be determined by the contact angle 9 at the wall. 
From the geometry we see that the polar angle (pi (> 0) at which curve (9) intersects 
the wall satisfies 

(13) 



a. 



r sin(0i + — ) = A sm — . 



Setting the angle between the interface and the wall (in three dimensions) to ^, we 
obtain, 

0. (14) 



e cos 6* + -— cos((/)i H — a) H — cos d{\ — x cos (p\)^ 
ad) 2 2 



Combining Eqs. (12)-(14) gives 
cos 6* 



C 



2 sin |q; 



a 



4 + 3sin^^- -(7r-20-a + sin2^)cot 



(15) 



Figure 3(a) shows cross-sectional profiles of the interface calculated using (12) (solid 
lines) for the case a = ^tt and 6 = j^^tt. The heights of the cross-sections are at 



zja = 1.5, 2, 2.5 and 3. The convergence to the circular profiles (shown by the 
dashed lines) is very fast as z/a increases. Figure 3(b) shows the vertical profiles of 
the interface for the same set of parameter values. Solid line gives the contact line 
of the liquid with the wall, while the dashed line is the profile on the bisector of the 
wedge. 

3. Dynamics 

When a thin tube with sharp grooves is made in contact with a liquid, the build- 
up of the equilibrium height consists of three stages: (i) an initial "rush" into the 
tube, which typically takes less than 10~^sec[3]. The flow in this period can be quite 
turbulent, (ii) viscous rising in the center part of the tube towards equilibrium, (iii) 
development of tongues. In the following we only consider stages (ii) and (iii). Here 
the dynamics is controlled by the viscous flow of the liquid which limits transport of 
matter needed to reach flnal equilibrium. For simplicity we assume that processes 
(ii) and (iii) are separated in time, though in practice there may well be an overlap. 
Our considerations follow closely previous work on wetting dynamics as reviewed by 
de Gennes[4] and by Leger and Joanny[5]. 

3.1. Viscous rise in a circular tube 

Let us flrst consider rising in a thin circular tube without grooves. The hydro- 
dynamic equation governing the viscous flow of the liquid in the tube is given by [6] 

-V{-+gz)+iyV^v = 0, (16) 

P 

where u is the viscosity coefficient of the liquid. The inertia term has been ignored. 
As usual, no-slip boundary condition on the wall is imposed. 

The driving force for the flow is the unbalanced pressure drop across the liquid-gas 
interface at height h, 

5p = pgh-4'jcos9/D = pg{h- heq), (17) 

where D is the diameter of the tube and heq = 2a^ cos 9/D is the equilibrium height 
of the interface[7]. The contact angle 9 has a weak dependence on velocity which we 
ignore here[4,5]. In the following we make the further approximation that v has only 
a vertical component. For an incompressible fluid, this implies that Vz depends only 
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on the horizontal coordinates {x^y). In this case Eq. (16) reduces to the Laplace 
equation 

{dl + dl)v, = 5p/{yph). (18) 

The solution to (18) is given by 

v. = \[6p/{uph)][r^-{D/2)\ (19) 

where r is the distance to the center of the tube. 

From Eq. (19) we obtain the flow rate per unit area 

, \ ,, I dxdy V, = --^(-)\ (20) 

Ti{D/2f J ^ %uph^2> ^ ' 

Identifying (20) with the velocity of the interface and using (17), we obtain, 

dh _ g .D.2heq^-h 

dt~8ry^2^ h ■ ^ ' 

This equation can be easily integrated to give 

^ Ml -1^) = ^ (-)'-, (22) 



/igq heq 8 COS 6* 2a to 

where to = ^/{9o) is a characteristic relaxation time of the liquid, typically less than 
10~^ sec. For h/heq <C 1 or t <C to{2a/D)^ the rising is diffusive, 

h = ^^(-f\-f'. (23) 

2cosi/2^^2a^ 4o^ 



It crosses over to an exponential decay 

h = h. 



"eq 



'-M-'-s^eO'B ''*' 



at t ~ to{2a/Df. 

The important result of the above calculation is that the typical rising time of 
the liquid is inversely proportional to the third power of the diameter of the tube. 
Another interesting phenomenon is that, although the equilibrium rise /leq is higher 
when the tube is thinner, the time it takes to reach a given height h <C /leq is actually 
longer, t{h) ^ toh'^ / (aD) . The diffusive regime exists only when the diameter is 
signiflcantly smaller than the capillary length. The approximation that v has only 



a vertical component is not correct close to the top of the liquid column or at the 
root of tongues. However, we expect that the qualitative behavior of viscous rising 
is captured by (22). 

3.2. Development of tongues 

In the above discussion we assumed that the equilibrium meniscus shape (but 
not the height) at the center of the tube is established at the end of stage (i) . The 
rising in stage (ii) is then governed by the unbalanced pressure drop, compensated 
by the viscous forces in the liquid. This description does not apply directly to the 
tongue region: the equilibrium interface extends to infinite height which can not be 
established in any finite time. A more plausible picture is that, at a given time t, 
there is a fully-developed tongue truncated at some height hmit). In addition, there 
can be an incipient tongue above hmit). 

A simple estimate for hm{t) can be made from Eq. (23). The equilibrium thick- 
ness of the tongue at hm is of the order of a? /hm [see Eq. (6)]. Assuming that the 
viscous flow up to this level is similar to the one in a tube of diameter D = o? /h^, 
we obtain from Eq. (23) 

/i^~a(t/to)'/'. (25) 

Equivalently, the time for the tongue to reach a height h is given by 

r ~ to{h/af. (26) 

It turns out that a more detailed calculation presented below gives the same result 
as the simple-minded estimate (26). Let us start with the approximate equation 

idl + dl)v, = <Piz), (27) 

where 

cl){z) = -{g + -d,p). (28) 

The term d'^Vz, which is omitted in (27), will be shown to be small. In addition to the 
no-slip boundary condition f ;2 = on the solid walls, we demand that the tangential 
stress vanishes on the liquid-gas interface, dnVz = 0. 

To flx the solution to (28), one has to specify the shape of the cross-sectional 
area. For simplicity we assume that all cross-sections of the tongue have the same 



shape, and are parametrized only by the hnear dimension R{z). In this case solution 
to (28) can be written as 

v,{x, y, z) = -v{x/R, y/R)R^(P, (29) 

where v is the solution to a dimensionless equation which vanishes on the walls and 
is positive inside the region of interest. We are now in a position to justify that d'^v^ 
is indeed negligible: it is smaller than (/> by a factor (R/z)'^, which is a small number 
in the tongue. 

The final step is to write down an equation for R{z, t), which can be done using 
the continuity equation for an incompressible liquid, 

dtA + d,U = 0. (30) 

Here A = AqR"^ is the cross-sectional area, U = f dxdy Vz{x,y) = —UqR'^c/) is the 
total velocity in the vertical direction, and Aq and Uq are positive numbers which 
depend only on the geometry. Combining (30) with (28) and using P — Pq — 7/-R, 
we obtain 

dtR^ = (t/oMo)c>4^i?4 ^ SLr'^q^rY (31) 

Equation (31) has a stationary solution (6). We now consider a linearized form 
of (31) around this solution. Writing R = (a^/2z)(l + 5), we obtain, 

t,dt6 = -dl5 - '^dj - ^6, (32) 

jO jO Ai 

where t\ = 8toAo/Uo. Equation (32) admits a scaling solution 
The scaling function A satisfies 

(PA (I 2 _l\(^_±^^_r. 

du^ 3 u du v? 

Thus the characteristic relaxation time for the build-up of the fully-developed tongue 
at height z is given by (26). 

In the long-time limit, 5 ^ 0, so that A(0) = 0. The second boundary condition 
may be chosen to be A(oo) = —1. A numerical integration of (34) then yields the 

9 



solution shown in Fig. 4. For tt ^ 1 we have A = —1 + 4u ^ + 0{u ^). In the 
opposite hmit w ^ 1 we have A ~ — 0.0317w^[l — j^u^ + 0{u^)]. 

4. Discussion and Conclusions 

In this paper we have shown that a liquid drop of finite extent is unstable against 
spreading in a groove whose opening angle is less than some critical value ac- Spread- 
ing to infinite heights takes place even when the liquid has to climb upwards against 
gravity. However, the thickness of the liquid tongue decreases with increasing height. 
In addition, the total amount of liquid in the tongue is finite: it is of the order of 
a^D where D is the diameter of the tube. 

Since the equilibrium thickness of the tongues is proportional to a? / z [Eq. (6)], 
at very high altitudes, microscopic interactions such as van der Waals and double- 
layer forces become important [1,3-5]. Taking a = 1 mm and the range of molecular 
interactions at 1/um, the macroscopic description used in this paper breaks down 
when the height exceeds 1 m. This regime is left for future investigation. 

In the analysis presented on the dynamics of the rising process, we assumed that 
the approach to equilibrium is controlled by viscous fiow of the liquid, which yields 
a t^/^-law for rising of the center meniscus and, at a later stage, a slower t^' ^-law 
for rising of the tongues. These results presumably apply to very clean surfaces 
where contact line hysteresis is negligible. For pure water at room temperature, the 
surface tension and kinematic viscocity coefficients are given by 7 ^ 73 dyn/cm and 
z/ ~ 10~^cm^/sec, respectively. This yields a = 0.38 cm and to = ^/{go) — 3 x 10~^ 
sec. Using Eq. (26), we find that the time it takes for the liquid tongue to reach 1 cm, 
10 cm, and 1 m are then 10~^ sec, 1 sec, and 10 min, respectively. In our experiments 
with colored solutions, rising seems to be much slower than that predicted above. 
The discrepancy at short times is expected from the fact that the Reynolds number 
R = vL/u is quite high there. At later times, contact-line pinning by dirts on the wall 
of the tube can infiuence the rising process dramatically. It is not clear, however, that 
the scaling laws should break down when the velocity of the contact line is sufficiently 
high. (For a discussion of contact line dynamics at low velocities see Refs. [4,5,10].) 

The analysis presented in this paper offers an explanation for the fast release of 
bubbles in tubes with sharp grooves. Such tubes may be more of a common occur- 
rence than rarity in living bodies, for which liquid transport is of crucial importance, 
and where capillary forces are of relevance. Our findings may also be of use in solving 
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engineering problems where trapping of bubbles is hazardous. We hope our prelim- 
inary study will inspire further theoretical and experimental investigations in this 
area, particularly on the dynamical side. 

We wish to acknowledge useful discussions with I. F. Lyusyutov and G. Nimtz. 
Upon completion of the work, we discovered that wetting in a wedge near the liquid- 
gas transition and in the absence of gravity has been discussed previously, and the 
importance of the critical groove angle otc has been realized in that context [8,9]. We 
thank H. Dobbs for bringing these work to our attention. The research is supported 
in part by the German Science Foundation (DFG) under SFB 341. 
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Figure Captions 

Fig. 1. Cross-sectional shape of tubes used in the experiment, (b)-(e) contain sharp 

grooves along the tube. The groove opening angles are: (b) a = ^tt; (c) 

a = 0; (d) a = ^tt; (e) a = iyr. 
Fig. 2. (a) A wedge of angle a formed by two vertical planes. Thin lines illustrate 

equilibrium capillary rise, (b) Horizontal cross-section of (a). Liquid is con- 

flned in the region OCBDO. 
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Fig. 3. Equilibrium profile of the tongue in an infinite wedge of opening angle a. = 
^TT. The contact angle is chosen at ^ = ^tt. (a) Horizontal cross-sections of 
the tongue (thin solid lines) at heights z/a = 1.5, 2, 2.5 and 3 (from right 
to left). Dashed lines give the zeroth order approximation, (b) The contact 
line (solid) on the wall and the vertical profile (dashed) on the bisector of the 
wedge. 

Fig. 4. The scaling function that appears in Eq. (33) 
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